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Abstract 



We investigate the singularity formation of a nonlinear nonlocal system. This non- 
' local system is a simplified one-dimensional system of the 3D model that was recently 

proposed by Hou and Lei in |16j for axisymmetric 3D incompressible Navier-Stokes 
equations with swirl. The main difference between the 3D model of Hou and Lei and 
the reformulated 3D Navier-Stokes equations is that the convection term is neglected in 
the 3D model. In the nonlocal system we consider in this paper, we replace the Riesz 
operator in the 3D model by the Hilbert transform. One of the main results of this 
paper is that we prove rigorously the finite time singularity formation of the nonlocal 
CN ■ system for a large class of smooth initial data with finite energy. We also prove the 

global regularity for a class of smooth initial data. Numerical results will be presented 
to demonstrate the asymptotically self-similar blow-up of the solution. The blowup rate 
of the self-similar singularity of the nonlocal system is similar to that of the 3D model. 

Key words: Finite time singularities, nonlinear nonlocal system, stabilizing effect of 
convection. 



1 Introduction 



The question of whether a solution of the 3D incompressible Navier-Stokes equations can 
develop a finite time singularity from smooth initial data with finite energy is one of the 
most outstanding mathematical open problems \20\ [23] . A main difficulty in obtaining 
the global regularity of the 3D Navier-Stokes equations is due to the presence of the vor- 
tex stretching term, which has a formal quadratic nonlinearity in vorticity. So far, most 
regularity analysis for the 3D Navier-Stokes equations uses energy estimates. Due to the 
incompressibility condition, the convection term does not contribute to the energy norm 
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of the velocity field or any LP (f < p < oo) norm of the vorticity field. In a recent paper 
by Hou and Lei [16J, the authors investigated the stabilizing effect of convection by con- 
structing a new 3D model for axisymmetric 3D incompressible Navier-Stokes equations with 
swirl. This model preserves almost all the properties of the full 3D Navier-Stokes equations 
except for the convection term which is neglected. If one adds the convection term back 
to the 3D model, one would recover the full Navier-Stokes equations. They also presented 
numerical evidence which supports that the 3D model may develop a potential finite time 
singularity. They further studied the mechanism that leads to these singular events in the 
3D model and how the convection term in the full Navier-Stokes equations destroys such a 
mechanism. 

In this paper, we propose a simplified nonlocal system for the 3D model proposed by 
Hou and Lei in |16| . The nonlocal system is derived by first reformulating the 3D model of 
Hou and Lei as the following two-by-two nonlinear and nonlocal system of partial differential 
equations: 

ut = 2uv + vAu, vt = (— A)~ 1 d zz u 2 + uAv, (1) 

where u = u 9 /r, v = ip z /r, and A = d z + d 2 + ^d r , and vP is the angular velocity 
component and ip® is the angular stream function respectively, r = \J x 1 + y 2 . By the 
partial regularity result for the 3D model [13], which is an analogue of the well-known 
Caffarelli-Kohn-Nirenberg partial regularity theory for the 3D incompressible Navier-Stokes 
equations [2], we know that the singularity can only occur along the symmetry axis, i.e. the 
z-axis. In order to study the potential singularity formation of the 3D model, it makes sense 
to construct a simplified one dimensional nonlocal system along the z-axis. One obvious 
choice is to replace the Riesz operator (— A)~ 1 d z by the Hilbert transform H along the z 
axis, and replace An by u zz , Av by v zz . This gives rise to our simplified nonlocal system: 

u t = 2uv + vu zz , v t = H(u 2 ) + uv zz , (2) 

where H is the Hilbert transform, 

(Hf) (x) = iRV. ^dy. (3) 

7T J-oo x-y 

In our analysis, we will focus on the inviscid version of the nonlocal system and relabel the 
variable 

u t = 2uv, v t = H(u 2 ), (4) 

with the initial condition 

u (t = 0) = u (x), v(t = 0) = v (x). (5) 

Note that the ID model ([2]) is designed to capture the dynamics of the 3D model (pQ) 
along the z-axis only. Thus, its inviscid model (j3|) does not enjoy the energy conservation 
property of the original model in the three-dimensional space. 

One of the main results of this paper is that we prove rigorously the finite time singularity 
formation of the nonlocal system for a large class of smooth initial data with finite energy. 
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As we will demonstrate in this paper, the blowup rate of the self-similar singularity of the 
nonlocal system ©-([5]) is qualitatively similar to that of the 3D model. The main result of 
this paper is summarized in the following theorem. 

Theorem 1.1 Assume that the support of uq is contained in (a,b) and that uq, vq € H 1 . 
Let 4>{x) = x — a and 



C = 4 / (f)(x)u Vq dx, /( 



oc 



dy 



vV+T 



If C > 0, then the solution of the nonlocal system (0)-([5[) must develop a finite time singu- 

AC ^ 1/3 
^3ir(b — a) 



larity in the H 1 norm no later than T* = ( - — — ) / ( 



oc • 



A similar result has been obtained for periodic initial data. 

The analysis of the finite time singularity for this nonlocal system is rather subtle. The 
main technical difficulty is that this is a two-by-two nonlinear nonlocal system. The key 
issue is under what condition the solution u has a strong alignment with the solution v 
dynamically. If u and v have a strong alignment for long enough time, then the right 
hand side of the u equation would develop a quadratic nonlinearity dynamically, which 
will lead to a finite time blowup. Note that v is coupled to u in a nonlinear and nonlocal 
fashion through the Hilbert transform. It is not clear whether u and v will develop such a 
nonlinear alignment dynamically. To establish such a nonlinear alignment, we need to use 
the following important property of the Hilbert transform: 

Proposition 1.1 Let 4> be a globally Lipschitz continuous function on R. For any f € 
W(R l ) n L 1 {R l ) and <j>f G L q {R l ) with 1 + 1 = 1, 1 < p, q< +oo, we have 

" <j>(x)f(x)Hf(x)dx = -L / + °° / + °° m ~_^ {y) f(x)f{y)dxdy. (6) 

3 ^ J — OO J — OO X y 

Using this property, we can identify an appropriate test function <p such that the time 
derivative of f u 2 (f)dz satisfies a nonlinear inequality. This inequality implies a finite time 
blowup of the nonlocal system. 

Proposition 11.11 should be a well-known property in the Harmonic Analysis literature. 
During the revision of our paper, we found that an identity which can be used to derive the 
special case (j> = x of Proposition 11.11 has been used in [10] , see also a recent paper [19] 0. 
However, we have not been able to find a proof for the general case stated in Proposition 
1.11 in the literature. For the sake of completeness, we provide a proof of Proposition 1 1 . 1 1 in 
Section 2. 

Another interesting result is that we prove the global regularity of our nonlocal system 
for a class of smooth initial data. Specifically, we prove the following theorem: 



1 We only learned about the work of [19] after the presentation of our work at the PIMS workshop on 
Hydrodynamics Regularity in August 2009. 
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Theorem 1.2 Assume that uq, vo £ H . Further we assume that uq has compact support 
in an interval of size 5 and vq satisfies the condition vq < —3 on this interval. Then the H 1 
norm of the solution of the nonlocal system remains bounded for all time as long as 

the following holds 

5 1/2 (\\vo x \\Li + l^\\u 0x \\h)< 1 -- (7) 



3 4 

Moreover, we have \\u\\l°° < Ce~ 3t , \\u\\gri < Ce~ 3t , and \\v\\hi < C for some constant C 
which depends on uq, Vq, and 5 only. 

In order to study the nature of the singularities, we have performed extensive numerical 
experiments for the nonlocal system with or without viscosity. Our numerical study shows 
that ||u||ioo(i) and |M|l°°(£) develop a finite time blowup with a blowup rate O Cjr-A, 
which is qualitatively similar to that of the 3D model [16J. Our numerical results also 
indicate that the solution of the inviscid nonlocal system seems to develop a one-parameter 
family self-similar finite time singularity of the type: 

u(x,t) = ^L.U^t), (8) 

v(x,t) = ^L^V&t), (9) 
f = x-x (t) 

? (T - i) 1 / 2 log(l/(T - t))V2 ' 1 > 

where xo(t) is the position at which u(x,t) achieves its maximum. The parameter that 
characterizes this self-similar blowup is the rescaled speed of propagation of the traveling 
wave defined as follows: 



Different initial data give different speeds of propagation of the singularity. One of the 
interesting findings of our numerical study is that by rescaling the self-similar variable £ by 
A , the different rescaled profiles corresponding to different initial conditions all collapse to 
the same universal profile. We offer some preliminary analysis to explain this phenomenon. 

Our numerical results also show that there is a significant overlap between the inner 
region of U and the inner region of V where V is positive. Such overlap persists dynamically 
and is responsible for producing a quadratic nonlinearity in the right hand side of the it- 
equation. The nonlinear interaction between u and v produces a traveling wave that moves 
to the righo Such phenomenon seems quite generic, and is qualitatively similar to that 
of the 3D model [16J. The only difference is that the 3D model produces traveling waves 
that move along the symmetry axis in both directions. It is still a mystery why the inviscid 
nonlocal system selects the scaling (jlOf) with the 1/2 exponent and a logarithmic correction. 
With the logarithmic correction, the viscous term can not dominate the nonlinear term 2uv 
in the equation. Indeed, when we add viscosity to the nonlocal system, we find that the 



2 If we change the plus sign in front of the Hilbert transform in the nonlocal system ([2]) to a minus sign, 
the nonlocal system would produce a traveling wave that moves to the left. 
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viscous solution still develops the same type self-similar finite time blowup as that of the 
inviscid nonlocal system. 

We remark that Hou, Shi and Wang |17] have recently made some important progress in 
proving the formation of finite time singularities of the original 3D model of Hou and Lei [16] 
for a class of smooth initial conditions with finite energy under some appropriate boundary 
conditions. The stabilizing effect of convection has been studied by Hou and Li in a recent 
paper |15] via a new ID model. Formation of singularities for various model equations for 
the 3D Euler equations or the surface quasi-geostrophic equation has been investigated by 
Constantin-Lax-Majda [6], Constantin [5], DeGregorio [BJ|9], Okamoto and Ohkitani |21| . 
Cordoba-Cordoba- Fontelos [7], Chae-Cordoba-Cordoba-Fontelos [4], and Li-Rodrigo |18j . 

The rest of the paper is organized as follows. In Section 2, we study some properties 
of the nonlocal system. In Section 3, we establish the local well-posedness of the nonlocal 
system. Section 4 is devoted to proving the finite time singularity formation of the inviscid 
nonlocal system for a large class of smooth initial data with finite energy. We prove the 
global regularity of the nonlocal system for a class of initial data in Section 5. Finally, 
we present several numerical results in Section 6 to study the nature of the finite time 
singularities for both the inviscid and viscous nonlocal systems. 

2 Properties of the nonlocal system 

In this section, we study some properties of the nonlocal system. First of all, we note that 
the nonlocal system has some interesting scaling property. Specifically, for any constants a 
and /3 satisfying a/3 > 0, the nonlocal system 

ut = auv, vt = f3Hu 2 (11) 

is equivalent to the system 

u t = 2uv, v t = Hu 2 (12) 
by introducing the following rescaling of the solution: 

u = u (x, jt) , v = fj,v (a?,7t) , (13) 

where 7 and fi are related to a and (3 through the following relationship: 

7 = W— , n = sgn[a)\ — . (14) 
V 2 V a 

Therefore, it is sufficient to consider the nonlocal system in the following form: 

u t = 2uv, vt = Hu 2 . (15) 

Moreover, if we replace the second equation vt = Hu 2 by v = Hu 2 and define w = u 2 , 
then our nonlocal system is reduced to the well-known Constantin-Lax-Majda model [6]: 

wt = AwHw. (16) 
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Before we end this section, we present the proof of Proposition 11.1 



Proof of Proposition ED Denote f e (x) = - f -^-dy, FJx) = 4>(x)f(x)]Jx) 

vr J\ x ~y\ >e x - y 

and f(x) = sup|/ £ (x)|. It follows from 1 
[3] that f e (x) — > Hf(x) a.e. x € R 1 and 



and f(x) = sup|/ e (x)|. It follows from the singular integral theory of Calderon-Zygmund 

e>0 



(18) 



LP < CpH/HiP. 

Therefore, we have F e (x) — > (j)(x)f(x)Hf(x) a.e. x € R 1 and |-F e (x)| < G(x), where 
G(x) = \4>(x)f(x)\f(x) satisfies 

\\G(x)\\ Ll < \\f(x)\\ LP \\cf>(x)f(x)\\ Lq 

< C p \\f(x)\\ LP U(x)f(x)\\ Lq <+oo. (17) 

Using the Lebesgue Dominated Convergence Theorem, we have 
(x)f(x)H(f)dx = lim / f(x)(j)(x)f e (x)dx 

= — lim / f{x)4>{x) [ t-^- dydx. 

7T e~>0j J\x-y\>e x ~ V 

Note that 

f | /W | ( I J d» < / i/wi ( / m^a.) a y 

J \J\x~y\>e \x-y\ J J \J e+\x-y\ ) 

< 2||0(x)/(x)|| L9 ||(6 + |x|)- 1 || iP I \f(y)\dy 

= C\\f(y)\\ L iU(x)f(x)\\ L9 <oo, 

for each fixed e > since / G L , 4>f G L q by our assumption, and C = ||(e + < oo 

for p > 1. Thus Fubini's Theorem implies that 

- / / I^ldydx=- [ [ f{ x )^ x )I^ldydx, (19) 



K J J\x-y\>eX~y IT J J \x-y\>e X ~V 

for each fixed e > 0. Furthermore, by renaming the variables in the integration, we can 
rewrite 1/2 of the integral on the right hand side of (|19p as follows: 

^11 f(x)f(y)^-dydx = -±- f [ f( x )f(y)^-dxdy, 

2^7 J\x-y\>e X ~V 27r J J\x-y\>e X ~V 

which implies that 

- / / mm^dydx = i- / / mm^^ldxdy. m 

* J J\x-v\>t x ~y 27r J J\x-y\>e x ~y 
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Since / € L l (R) and <j>{x) is globally Lipschitz continuous on R, it is easy to show that 

x-y 

Using the Lebesgue Dominated Convergence Theorem, we have 

h hS L., /w/m^o^ =± fj mm *J^m Axi , (21) 

Proposition 11.11 now follows from (|18p -(|2i p . □ 

We remark that Proposition 11.11 is also valid for periodic functions. Recall that for 
periodic functions (with period 2tt) the Hilbert transform takes the form: 

(Hf) (x) = i-P.V. jT f (y) cot (^*\ dy. (22) 

For the sake of completeness, we state the corresponding result for periodic functions below: 

Proposition 2.1 Let <p be a periodic Lipschitz continuous function with period 2tt. For 
any periodic function f with period 2tt satisfying f € L p ([0,27r]) and 4>f £ L 9 ([0,27r]) with 
| + | = 1 ; 1 < p, q < +oo, we have 

2T 1 f 27T (** tJ ., x j., (*-V 



<j>{x)f{x)Hf{x)6x = — J J (<P(x) - 0(y)) cot y-^) f(x)f(y)dxdy. (23) 

The proof of Proposition 12.11 goes exactly the same as the non-periodic case. We omit 
the proof Here. 

Remark 2.1 As we see in the proof of Proposition HOI the key is to use the oddness of the 
kernel in the Hilbert transform. The same observation is still valid here: 

I I f(x)f(y)(j)(x) cot [ ^— — ) dydx 



27r ■/ -/[O^tt] 2 , \x-y\>e ' V - 

f(x)f(y)<j>(y) cot — — ) dxdy, 



2tt 



[0,2tt] 2 , \x-y\>e 



by renaming the variables in the integration. 



3 Local well-posedness in H 

In this section, we will establish the local well-posedness in Sobolev space H l . 

Theorem 3.1 (Local well-posedness) For any uo,vq € H , there exists a finite time T = 
r(||«o||i?ij ll^olli^ 1 ) > such that the nonlocal system (0)-([5|) has a unique smooth solution, 
u,v £ C 1 ({0,T); H 1 ^) for < t < T. Moreover, if T is the first time at which the solution 
of the nonlocal system ceases to be regular in H 1 and T < oo, then the solution must satisfy 
the following condition: 

T 

(IM|l°° + di = +oo. (24) 
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Remark 3.1 We remark that the condition \2J$ is an analogue of the well-known Beale- 
Kato-Majda blowup criteria for the 3D incompressible Euler equation JT^. 

Proof To show local well-posedness, we write the system as an ODE in the Banach space 
X := H 1 x H 1 : 

U t = F(U), (25) 

where U = (u,v), F(U) = (2uv , H(u 2 )). As H (R) is an algebra, F maps any open set 
in X into X, and furthermore F is locally Lipschitz on X. Local well-posedness of (jH)-© 
then follows from the standard abstract ODE theory such as Theorem 4.1 in |20j. 

The blow-up criterion ()24|) follows from the following a priori estimates. Multiplying 
the u-equation by u and the f-equation by v, and integrating over R, we obtain 

— f u 2 dx = 4 f u 2 vdx < 4||u||i,oo f u 2 dx, (26) 



dt 



and 



d_ 

dt 



J v 2 dx = 2 j vHu 2 dx = -2 j (Hv)u 2 dx<2\\u\\ L oo j \Hv\udx 

< \\u\\ L °° I / u 2 dx+ / v 2 dx ) . (27) 



Similarly, we can derive L 2 estimates for u x and v x as follows: 
di 



u x dx = 4 / (vu x + uv x u x )dx 

< m\^ju 2 x dx + 2 M ^j\u 2 x + v 2 x )dx, (28) 



and 



d_ 

dt 



v 2 dx = 4 / v x H (uu x ) dx 
= 4 J (Hv x )uu x dx 



< 2\\u\\ L <*> J (v£ + v£)dx. (29) 
Summing up the above estimates gives 

^ (IMIffi + M\hi) < C(\\u\\v*> + \\v\\ L ™) (\\u\\ 2 H i + Hull^i) . (30) 
We see that the regularity is controlled by the quantity 

IMIz,°° + IMIz,°° • 

If / T (ll u lli°° + II^IIl 00 )^ < °°) then it follows from ([30]) that ||n||jji + must remain 

finite up to T. Therefore, if T is the first time at which the solution blows up in the 
i? 1 -norm, we must have 

r-T 

(IMIl°° + IM|l°°) dt = +oo. (32) 

□ 
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4 Blow up of the nonlocal system 



In this section, we will prove the main result of this paper, that is the solution of the 
nonlocal system will develop a finite time singularity for a class of smooth initial conditions 
with finite energy. We will prove the finite time singularity of the nonlocal system as an 
initial value problem in the whole space and in a periodic domain. 



4.1 Initial Data with Compact Support 

We first consider the initial value problem in the whole space and prove the finite time blow 
up of the solution of the nonlocal system (H])-© for a large class of initial data -uo that have 
compact support. 

For the sake of completeness, we will restate the main result below: 

Theorem 4.1 Assume that the support of uq is contained in (a,b) and that uo, vq € H 1 . 
Let 4>(x) = x — a and 



C = 4 / 4>(x)uq vo dx, 



dy 



Vy 3 + i' 



If C > 0, then the solution of the nonlocal system (^-([5]) must develop a finite time singu- 

( AC \ -1 / 3 

larity in the H 1 norm no later than T* = I — — I 1^. 

\3ir(b — a) 2 J 



Proof By Theorem l3,H we know that there exists a finite time T = T (\\uq\\jji, H^oH-H 1 ) > 
such that the nonlocal system ([I])-© has a unique smooth solution, u,v G C 1 ([0, r);^ 1 ) 
for < t < T. Let T* be the largest time such that the nonlocal system with initial 
condition uq and vq has a smooth solution in H . We claim that T* < oo. We prove this 
by contradiction. 

Suppose that T* = oo, i.e. that the nonlocal system has a globally smooth solution in 
H 1 for the given initial condition no and vo. Using (jU, we obtain 

{u 2 ) tt = A(u 2 v)t = 8u t uv + 4u 2 v t = 4{u t ) 2 + Au 2 H{u 2 ). (33) 

Multiplying 4>(x) to both sides of the above equation and integrating over [a,b], we have 
the following estimate: 



dt 2 



f'b rb 

4>(x)u 2 (x,t)dx = 4 / (j)(x)(u t ) 2 dx + 4 / (f)(x)u 2 H (u 2 )dx 

J a J a 

/■b 

> 4 / 4>(x)u 2 H(u 2 )dx . (34) 



Note that the support of u(x,t) is the same as that of the initial value uq. Proposition 11.1 
implies that 

b poo 

(j)(x)u 2 H(u 2 )dx = / <j)(x)u 2 H(u 2 )dx 
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2vr 

Combining (|34l) with (i35j) . we get 



I / fb \ 2 

u 2 (x,t)dx) . (35) 



dt 2 



rb \ 2 

(x)u 2 (x,f)dx > - ( / u 2 (x,£)dx . (36) 
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As we can see, Proposition 11,11 plays an essential role in obtaining the above inequality, which 
is the key estimate in our analysis of the finite time singularity of the nonlocal system. 
By the definition of (f>, we have the following inequality: 

b ^ rb 

u 2 (x,t)dx > / (f)(x)u 2 (x, t)dx. (37) 

b-a J a 

Combining (|36|) with (|37p . we obtain the following key estimate: 

^ J ^(x)u 2 (x, t)dx > ^ b 2 _ a)2 (J </>(x)v 2 (x, t)dx\ (38) 

Denoting F(t) = (fi(x)u 2 (x , t)dx we obtain the ODE inequality system 

2 



Ftt > 



ir(b-a) 2 



F 2 , F t (0) = C > 0, F(0) = I 0ug > 0. (39) 

J a 



Since F t (0) = C > 0, integrating (jpj) from to t gives F t (t) >0 for all t > 0. Denote 
F(t) = F(f) -F(0). Then we have F(t) > for t > 0, F t > and F(0) = 0. Since F(0) > 
and F(t) > 0, it is easy to show that F satisfies the same differential inequality (]39p as F. 
Therefore we can set F(0) = in the following analysis without loss of generality. 
Multiplying Ft to F# > ^z^p an d integrating in time, we obtain 



dF ^ 4 
"eft - V 3vr(6 - a 



> A / Q _„, ^ 2 F3 + C 2 . (40) 



It is easy to see from the above inequality that F must blow up in a finite time. Define 

x dy (M{b-a) 2 C 2 ^ 1/3 



I(x) = / -f=, J 
Jo vV + 1 



Integrating (|40p in time gives 
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Observe that both / and F are strictly increasing functions, and I{x) is bounded for all 
x > while the right hand side of (|4ip increases linearly in time. It follows from (|4ip that 
F{t) must blow up no later than 

(42) 

This contradicts with the assumption that the nonlocal system has a globally smooth solu- 
tion for the given initial condition no and vq. This contradiction implies that the solution of 
the nonlocal system (|4])-([5]) must develop a finite time singularity in the H 1 norm no later 
than T* given by (I42|) . This completes our proof of Theorem 14.11 □ 



4.2 Periodic Initial Data 

In this subsection, we will extend the analysis of finite time singularity formation of the 
nonlocal system to periodic initial data. Below we state our main result: 

Theorem 4.2 We assume that the initial values uq,Vq are periodic functions with period 
2ir and the support of uq is contained in (a,b) C (0, 2ir) with b — a < ir. Moreover, we 
assume that uq, vq G H [0, 2ir]. Let 4>{x) be a 2ir-periodic Lipschitz continuous function 
with 4>(x) = x — a on [a,b], and 



C = 4 4>(x)ul V dx, loo = 

J a JO 



dy 



If C > 0, then the solution of the nonlocal system ^-(j^) must develop a finite time singu- 

/ 4Ccos(^)V 1/3 



larity in the H 1 norm no later than T* 



3tt(6 - af 



Proof As in the proof of Theorem l4.ll we also prove this theorem by contradiction. Assume 
that the nonlocal system with the given initial condition uq and vq has a globally smooth 
solution in H . As before, by differentiating (j3J) with respect to t, we obtain the following 
equation: 

(u 2 ) tt = 4( Ut ) 2 + 4 U 2 H(u 2 ). (43) 

Multiplying cj>(x) to both sides of the above equation, integrating over [0, 2tt] and using 
Proposition 12.11 we obtain the following estimate: 

d 2 f b d 2 f 2lT 

—t; \ 4>(x)u 2 (x,t)dx = / (f>(x)u 2 (x, t)dx 

dt z J a dt 2 J 

n27r r 2n 



= 41 (j)(x)(u t ) 2 dx + 4 / (j){x)u 2 H{u 2 )dx 
Jo Jo 

[■2lT 

> 4 <f)(x)u 2 H(u 2 )dx 
Jo 
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*Jo I u2 ( x '*) n2 ^'*)^( a; ) -<^)) cot (^y^) d y dx 



a J a 



u 2 (x, t)u 2 (y, t)(x — y) cot ( — - — I dydx 



> — (^J u 2 (x,t)dx^j 



2 



(44) 



where M = min xcot(x/2). Since b — a < tt, we have 

— (6— a)<x<b— a 



xcos(x/2) . , , / b — a\ . 
M= min / / / > min 2 cos(x/2) = 2 cos > 0. (45) 

-(6-o)<a;<6-a sin(a;/2) -{b-a)<x<b-a \ 2 J 

Now, following the same procedure as in the proof of Theorem 14.11 we conclude that the 
solution must blow up no later than 



1/3 / 4Ccos^ ^' 1/3 
V6^(6-a) 2 J Ioo ~ ^3^(6-0) 

This contradicts with the assumption that the nonlocal system with the given initial condi- 
tion uq and vq has a globally smooth solution. This contradiction implies that the solution 
of the nonlocal system (H])-© must develop a finite time singularity in the H 1 norm no 
later than T* given by (|46p . This completes the proof of Theorem 14.21 □ 

Remark 4.1 We can also prove the finite time blowup of a variant of our nonlocal system 

u t = 2uv, v t = -H(u 2 ), (47) 

by choosing the test function 4>(x) = b — x. It is interesting to note that while the solution 
of produces traveling waves that propagate to the right, the solution of fl^7| ) produces 
traveling waves that propagate to the left. 

Remark 4.2 Our singularity analysis can be generalized to give another proof of finite time 
singularity formation of the Constantin-Lax-Majda model without using the exact integra- 
bility of the model. More precisely, we consider the Constantin-Lax-Majda model: 

u t = uH(u), 
u(t = 0) = uq(x), 



By choosing 4>{x) = x — a and following the same procedure as in the proof of Theorem \4-l\ 
we can show that if uq is smooth and has compact support, supp uq = [a,b] and uo(x) > 
on (a,b), then the L\ norm of the solution of l[48\ ) must blows up no later than 

T* = f y"-°> 2 . (49) 
<P(x)uq dx 
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Below we will give a different and simpler proof of the finite time blowup for the Constantin- 
Lax-Majda model. 

Multiplying (j)(x) to both sides of equation |^<§| ), integrating over the support (a,b), and 
using Proposition we obtain 



1 ' ^ 2 



(x — a)udx = (x — a)uH{u)dx = — / udx . (50) 

dt J a J a 27r \Ja J 

As J^(x — a)udx < (b — a) udx due to u > for x € [a, b], setting F(t) = J (x — a)udx 
we have 

1 rb 

?2 



F * ^ X2 F ' = / H*)u dx > 0. (51) 

2vr(6 - a) 2 J a 

This leads to 

which implies the finite-time blowup of F no later than T* = ? n ^ b a ^ — . 

J a <t>(x)u dx 

Similar result can be obtained for periodic initial data following the same analysis of 
Theorem \4-£\ 



5 Global regularity for a special class of initial data 

In this section, we will prove the global regularity of the solution of our nonlocal system for 
a special class of initial data. Below we state our main result in this section. 

Theorem 5.1 Assume that uq, vq £ H . Further we assume that uq has compact support 
in an interval of size 5 and satisfies the condition vq < —3 on this interval. Then the H 1 
norm of the solution of the nonlocal system ^-(^ remains bounded for all time as long as 
the following holds 

^(iboxlUa + ^IK.Ili,) <~ (53) 

Moreover, we have \\u\\l°° < Ce~ 3t , \\u\\hi < Ce -3 *, and \\v\\hi < C for some constant C 
which depends on uq, vq, and 5 only. 

Proof Note that (J53]) implies that e) 1//2 ||t;o x || L 2 < \ which gives -4 + 25 1 / 2 ||v 0a ; Hz, 2 < 
—3.5. By using an argument similar to the local well-posedness analysis, we can show that 
there exists Tq > such that ||u||#i and ||u||#i are bounded, v < — 2 on supp(u), and 
25 1 / 2 \\v x \\ L 2 < 1 for < t < T . 

Let [0,T) be the largest time interval on which ||u||z/i and are bounded, and both 

of the following inequalities hold: 

v < — 2 on supp(n) and 2<5 1//2 \\v x \\ L 2 < 1. (54) 

We will show that T = oo. 
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We have for < t < T that 



d_ 

df 



/»> = 4/ (ra l + „„ I „ I )d,<-8/«^ + 4||„|| 1 „|K|| IS ||„ I || 1 , (55) 



Observe that supp(u) = supp(uo) for all times. Let £1 = supp(uo). Since supp(u) has length 
5, we can use the Poincare inequality to get 



\U\\L° 



<S 1/2 \WxhHn) = 5 1 / 2 \\u x \\ L 2. (56) 



Therefore we obtain the following estimate: 
d 



iLU u x \\ L 2 < -A\\u x \\ L 2 + 2S 1 ' 2 \\V X \\ L 2 \\u x \\ L 2 
at 



^ -4 + 25 1/2 \\v x \\ L 2) \\u x \\ L 2 < -3\\u x \\ L 2. (57) 
Thus we have for < t < T that 

\\ux\\l 2 < \\ u 0x\\L 2e ~ 3t ■ ( 58 ) 

On the other hand, we have that 

J v 2 x dx = A J v x H (uu x )dx < A\\v x \\ L 2\\uu x \\ L 2 

< 4\\u\\ L °°\\v x \\ L 2\\u x \\ L 2 < 45 1/2 ||t> a; || L 2||u x .|| 2 2 , (59) 



d 

dt 



where we have used the property that ||-f^(/)||i2 < \\f\\i 2 and the Poincare inequality 
Now using ([58]) . we get 

^KHl2 < 25 1 / 2 \\u x \\ 2 L2 < 2<5 1 / 2 ||n 0:c |||2e- 6 '. (60) 
As a consequence, we obtain for < t < T that 

IMIl2 < \\vox\\l* + ^ 1/2 \\uq x \\ 2 L 2. (61) 

Now observe that at the left end of the support of u, vt = Hu 2 is always negative. Since 
vo < —3 on the support of u, we conclude that v(x,t) < —3 at the left end of the support 
of u for all times. Now, we apply the Poincare inequality in the support of u and use (|6ip 
to obtain 



v < -3 + ^/211^11^^ < s + 5 1 / 2 (\\v 0x \\l* + ls l/2 \\u x\\ 2 L 2^ , 

on supp(tt) for all t G [0,T). 

Next, we perform L 2 estimates. We can easily show by using vt = Hu 2 that 



(62) 



~2dt 



J v 2 dx = j vHu 2 dx < \\v\\ L 2 \\u 2 \\ L 2 < \\v\\ L 2 \\u\\ 2 L acS 1 ^ 2 , 
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which gives 

^ II II ^ c-l /2 1 1 ||2 

^||w|| L 2 < 6 i \\u\\ Lao . 

It follows from (|56|) and (EHJ) that 

\\u\\ L °o < <5 1/2 ||u 0a; || L 2e- 3 *. (63) 



Therefore, we obtain 



which implies 



dt 1 



v\\& < S\\u 0x \\ 2 L2 e 6 *, 



v \\l2 < \\vo\\l* + l$\\ u 0x\\ 2 L 2, (64) 



6 

for < t < T. 

Similarly, using v < —2 on the support of u, we can easily show that 

IMIl 2 < IKHi^e -4 *, (65) 

for < t < T. 

To summarize, we have shown that and ||i>||_h-i are uniformly bounded for < 

t <T, and 

\\v x \\l 2 < HuOallia + — 5 1/2 ||^0a; 11^2 , (66) 

and 

v<-Z + 5 l l 2 \\v x \\ L 2 < -3 + 5 1 / 2 (\\v 0x \\ L 2 + U l l 2 \\u Qx \\ 2 L ^j , (67) 

on supp(it) for < t < T. 

By our assumption on the initial data, we have 

s 1/2 (\M\l* + V /2 IK,||| 2 ) <\. (68) 

Therefore, we have proved that if 

v < -2 on supp(u) and 25 1 / 2 ||i? a; || L 2 < 1, (69) 
< t < T, then we actually have 

v < —2.75 on supp(-u) and 25 l/2 \\v x \\ L 2 < 0.5, (70) 



< t < T. This implies that we can extend the time interval beyond [0, T) so that (|69p 
is still valid. This contradicts the assumption that [0, T) is the largest time interval on 
which (|69p is valid. This contradiction shows that T can not be a finite number, i.e. (|69p 
is true for all times. This in turn implies that and ||i>||#i are bounded for all times. 

Moreover, we have shown that both and decay exponentially fast in time and 

IMI/fi is bounded uniformly for all times (see ([63]) . (|58|) . (fBTj) . (jMj) and ([65])). This proves 
Theorem 15.11 □ 
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6 Numerical Results 



In this section, we perform extensive numerical experiments to study the nature of the 
singularities of the nonlocal system. Our numerical results demonstrate convincingly that 
the nonlocal system develops asymptotically self-similar singularities in a finite time for 
both the inviscid and the viscous nonlocal systems. 



6.1 Set-up of the Problem 

In our numerical study, we use the followingnonlocal system without the factor of 2 in front 
of the nonlinear term uv in the u-equation □: 

ut = uv + vu xx , (71) 
v t = H{u 2 ) + vv xx . (72) 

We study the nonlocal system for two types of initial data. The first type of initial data has 
compact support. The second type of initial data is periodic. The nature of the singularities 
for these two types of initial data is very similar. In the case of periodic data, we can use 
FFT to compute the Hilbert transform. This enables us to perform our computations with 
a very high space resolution. 

Below we describe the initial data that we will use in our numerical experiments. We 
choose three different initial conditions. The first initial condition has compact support 
which lies in Vt = [0.45, 0.55] and vq = 0. Within the compact support f2, uq is given by 

Initial Condition I: u$ = exp | 1 — ( 1 — ( - ^ ) ] | for vq = 0. 

We call this Initial Condition I. The largest resolution we use for Initial Condition I is 
N = 16,384. The timestep is chosen to be At = 10 — 3 /||ti||i°o in order to resolve the 
maximum growth of 

The last two initial conditions are periodic with period one. They are given as follows: 

Initial Condition II: uq = 2 + sin(27rx) + cos(47rx), v$ = 0, 

Initial Condition III: uq = ; -, vq = 0. 

1.2 + cos(2"7T2;) 

We call them Initial Condition II and Initial Condition III respectively. The largest reso- 
lution that we use for these two periodic initial conditions is N = 262, 144 = 2 18 , and the 
timestep is chosen to be At = 10~ 3 /||u||l°o . 

We use the fourth order classical Runge-Kutta method to discretize the inviscid nonlocal 
system in time. For the viscous nonlocal system, we consider only periodic initial data since 
the solution will not have compact support anymore. In order to remove the stiffness of the 



3 As we have shown in Section 2, dropping this factor only changes the scaling of the solution. 
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time discretization due to the viscous term, we first apply Fourier transform to the nonlocal 
system to obtain 



d t u(k, t) = (uv) (k,t)- vk 2 u{k, t) , (73) 
d t v(k,t) = -isgn(k)(u?)(k,t) -vk 2 v{k,t), (74) 

where u(k,t) is the Fourier transform of u and k is the wave number. We then reformulate 
the viscous term as an integral factor 

^(e vfca *«(M)) = e^M(M), (75) 

^(e"* a *«(M)) = -isMky kH &)(k,t). (76) 

Now we apply the classical Runge-Kutta method to discretize the above system in time. 
The resulting time discretization method will be free of the stiffness induced by the viscous 
term. 

For periodic initial data, we use the spectral method to discretize the Hilbert transform 
by using the explicit formula H(k) = — isgn(fc). For initial data of compact support, we use 
the well-known alternating trapezoidal rule to discretize the Hilbert transform which gives 
spectral accuracy. For the sake of completeness, we describe the method below, see also 
[22]. Let Xj = jh the grid point and h > is the grid size. The alternating trapezoidal rule 
discretization of the Hilbert transform is given by the following quadrature: 

H{f)( Xi )= £ — — 2h - ( 77 ) 

(j—i)odd 

Therefore, our numerical method has spectral accuracy in space and and fourth order ac- 
curacy in time. The high order accuracy of the method plus high space resolution and 
adaptive time-stepping is essential for us to resolve the asymptotically self-similar singular 
solution structure of the nonlocal system. 



6.2 Asymptotically self-similar blowup of the inviscid nonlocal system 

In the singularity analysis, we have proved that the nonlocal system must develop a finite 
time singularity for a large class of initial data. However, the singularity analysis does not 
tell us the nature of the singularity. Understanding the nature of the singularity is the main 
focus of our numerical study. Our numerical results show that for all three initial conditions 
we consider here, they develop asymptotically self-similar singularities in a finite time. The 
numerical evidence of self-similar singularities is quite convincing for all three initial data 
that we consider. As is the case for the original 3D model, the mechanism of forming such 
self-similar blowup of the nonlocal system is due to the fact that we neglect the convection 
term in our model. As it is demonstrated in [15\ [T6] , the convection term tends to destroy 
the mechanism for generating the finite time blowup in the ID or 3D model. Indeed, a 
recent numerical study shows that the 3D incompressible Euler equation does not seem to 
grow faster than double exponential in time \13\ 112] , 
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We use the following asymptotic singularity form fit to predict the singularity time and 
the blowup rate: 

NU- = (t^> ( 78 ) 

where T is the blowup time. We find that near the singularity time, the inverse of 
is almost a perfect linear function of time, see Figure [TJ To obtain a good estimate for the 
singularity time, we perform a least square fit for the inverse of ||u||l°o. We find that a = 1 
gives the best fit. The same least square fit also determines the potential singularity time 
T and the constant C. 

To confirm that the above procedure indeed gives a good fit for the potential singularity, 
we plot || U l l f-yQ CIS £l function of time with a sequence of increasing resolutions against the 
asymptotically form fit for the three initial conditions we consider here. In Figure [H we per- 
form such comparison for Initial Condition I with a sequence of increasing resolutions from 
N = 4096 to N = 16384. We can see that the agreement between the computed solutions 
and the asymptotically fitted solution is excellent as the time approaches the potential sin- 
gularity time. In the lower box of Figure [TJ we plot Hujloo computed by our adaptive method 
against the form fit C/(T-t) with T = 2.36752830915169 and C = 1.67396437016231. The 
computed solutions and the asymptotically fitted solution are almost indistinguishable. This 
asymptotic blowup rate is qualitatively similar to that of the 3D model |16j . 

We have also performed a similar comparison between the computed ||it||i,°° and the 
asymptotically fitted solution for Initial Conditions II and III in Figures [2] and Figures [3] 
respectively. For these two periodic initial conditions, we can afford even higher resolutions 
ranging from N = 2 14 to N = 2 18 . Again, we observe excellent agreement between the 
computed solutions and the asymptotically fitted singular solution. 

After we obtain an estimate for the singularity time, we can use it to look for a dynam- 
ically rescaled profile U(£, t), V(£, t) near the singularity of the form 

u(x, t) = — i— U ( X ~ Xo( ^ , t) , ast^T, (79) 



T-t V ( T ~ if ' 
ttfat) = -J—V ( X ~ X °^ ,t), ast^T, (80) 



t \{T- tf 

where T is the predicted blowup time in the singularity form fit (|78p . (3 is a parameter to 
be determined, and xo(t) is the location in which \u\ achieves its global maximum at t. 

Again, we use a least square fit to determine j3. Our numerical study indicates that 
/3 = 4 with a logarithmic correction. More precisely, we find that the dynamically rescaled 
variable £ has the form: 

c = x ~ x o( f ) / 8n 

? (T-ty/ 2 iog{i/(r-t)y/ 2 ' { } 

In terms of this rescaling variable £, we define the dynamically rescaled profiles U(£,t) and 
V(£, t) through the following relationship: 

u(x,t) = J—Ufot), (82) 
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Figure 1: Top: The inverse of (black) versus the asymptotic fit (red) for Initial 

Condition I with v = 0. The fitted blowup time is T = 2.36752830915169 and the scaling 
constant is C = 1.67396437016231. Bottom: ||u||t<x> (black) versus the asymptotic fit (red) 
for Initial Condition I. 



19 



,x 10" 



N=2" 




0.78 0.7801 0.7802 0.7803 0.7804 0.7805 0.7806 0.7807 0.7808 0.7809 

Time 



Figure 2: The inverse of ||ii||z,°° (black) versus the asymptotic fit (red) for Initial Condition 
II with v = 0. The fitted blowup time is T = 0.780894805082166 and the scaling constant 
is C = 1.68253514799506. 
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Figure 3: The inverse of ||it||z,°° (black) versus the asymptotic fit (red) for Initial Condition 
III with v = 0. The fitted blowup time is T = 0.569719056780405 and the scaling constant 
is C = 1.68293676812485. 
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-1 

-0.1 -0.05 0.05 0.1 

Figure 4: Rescaled profiles U and V for Initial Condition I with v = at three different 
times: t = 2.36710445318745, 2.36743324526419 and 2.36750705502071, the corresponding 
maximum values ofu are 3948, 17617 and 78422 respectively. Blue: profile of u; Red: profile 
of v . 

v(x,t) = ^L-Vitt). (83) 

In Figure 01 we plot the self-similar profiles U and V at three different times for Initial 
Condition I. We can see that the rescaled profiles for these three different times agree with 
one another very well. From Figure HI we can see that there is a significant overlap between 
the inner region of U and the inner region of V where V is positive. Such overlap persists 
dynamically and is responsible for producing a quadratic nonlinearity in the right hand side 
of the u-equation, which has the form 2uv. On the other hand, we observe that the position 
at which u achieves its global maximum is not in phase with the position at which v achieves 
its global maximum. In fact, the positive part of V always moves ahead of U. This is a 
consequence of the property of the Hilbert transform. As a result, the nonlinear interaction 
between u and v produces a traveling wave that moves to the right. Such phenomena seem 
quite generic. We observe the same phenomena for all three initial conditions for both the 
inviscid and the viscous models. This phenomenon is also qualitatively similar to that of 
the 3D model |16j . 

The strong alignment between the rescaled profile of u and v is the main mechanism 
for the solution of the nonlocal system to develop an asymptotically self-similar singularity 
in the form given by (|8ip and (|82p -(|83p. We observe essentially the same phenomena for 
Initial Conditions II and III, see Figure [5j 

It is interesting to see how the different rescaled profiles corresponding to different initial 
conditions are related to one another. In Figure [5] (top), we put three profiles from three 
different initial conditions together. The profile from Initial Condition III is the widest while 
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the profile from Initial Condition II is narrower than that from Initial Condition III. The 
profile from Initial Condition I is the narrowest of the three initial conditions. But what is 
amazing is that they can match each other very well by rescaling the variable variable £. To 
match the three rescaled profiles, we keep the profile from Initial Condition III unchanged. 
In order to match the profile from Initial Condition III, we change the profile from Initial 
Condition II by rescaling £ — >• £/1.58, and change the profile from Initial Condition I by 
rescaling £ — > £/19.5. As we can see from Figured] (bottom), the three rescaled profiles is 
almost indistinguishable. 

To gain some insight into this phenomenon, we perform some analysis of the self-similar 
solutions. We assume that the self-similar profiles converge to a steady state as t — > T. 

u(x, t) — > 7p-—j: U (£; A) , as t -> T (84) 
v(x,t) -> j^V&X), ast^T, (85) 

where A = lim f(T - t ) 1/2 ^ t x ^ t )\ • 

If we neglect the logarithmic correction in £ and substitute above equations into the 
nonlocal system, we obtain the equations for U and V as follows: 

U + ^U^-XU^ = UV, (86) 
V + PgVt-Wt = H(U 2 ). (87) 

Let Ui(£),Vi(£) be the solution of the self-similar system ([86]) . ([87]) corresponding to 
A = 1, then the solution for A / 1 can be obtained by using the following rescaling of the 
self-similar variable £: 

17(£; A) = tMA- 1 ?), (88) 
y(£;A) = ^(A-^). (89) 

The profiles that are obtained from different initial conditions have different A, but they 
can match each other by rescaling £. This may explain why we can match different rescaled 
profiles corresponding to different initial conditions by rescaling £. 

6.3 Asymptotically self-similar blowup of the viscous nonlocal system 

In this subsection, we perform computations to investigate the finite time singularity of 
the viscous nonlocal system. In our computations, we choose the viscosity coefficient to 
be v = 0.001. Notice that the solution of the viscous nonlocal system can not keep the 
compact support, so we only perform our numerical study for Initial Conditions II and III 
which are periodic. The computational settings are the same as those in the inviscid case. 
We use the same asymptotic singularity form fit as in the inviscid model, i.e. 

C 

\W\\L°° = {T _ t)a , (90) 
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Figure 5: The self-similar profiles for Initial Conditions I-III respectively (y = 0). Top: 
The original profiles for u for Initial Conditions I-III; Bottom: The rescaled profiles. Black: 
Initial Condition I; Red: Initial Condition II; Blue: Initial Condition III. 
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where T is the blowup time. In Figure [6] and Figure [71 we plot ||u||7jo versus the asymptotic 
singularity fit. We can see that as we increase resolutions from N = 2 14 to N = 2 18 , ||u||7oo 
converges to the asymptotic fit which is almost a perfect straight line. This suggests that 
a = 1. From these numerical results, we can see that adding viscosity with u = 0.001 does 
not prevent the solution from blowing up and does not change the qualitative nature of the 
singular solution, although it postpones the blowup time. 

Next, we study the rescaled profiles of the asymptotically self-similar solutions of the 
viscous nonlocal system. We look for a dynamically rescaled profile £/(£, t), V(£, t) near the 
singularity of the form 

u(x, t) = — !— U ( X ~ X °® , t] , asi^T, (91) 



T-t \(T-t)P 
v ( x , t ) = — L_ v( X ~ Xo ® ,t), ast^T, (92) 



T-t \(T -t 

where T is the predicted blowup time in the singularity form fit (j90|) , /3 is a parameter to be 
determined, and Xo(t) is the location in which \u\ achieves its global maximum at t. Again, 
we use a least square fit to determine f3 and find that (3 = 1/2 with a logarithmic correction. 
In Figures [8] and [9j we plot the rescaled profiles of the asymptotically self-similar solution 
for Initial Conditions II and III respectively. The dynamically rescaled variable £ has the 
same form as that of the inviscid nonlocal system, i.e. 

= x - X Q (t) 

? (T-ty/ 2 \o g (i/(T-t)y/ r { 1 

In Figure [8j we plot the self-similar profiles U and V at three different times for Initial 
Condition II. We can see that the rescaled profiles for these three different times agree with 
one another very well. As in the inviscid case, we observe that there is a significant overlap 
between the inner region of U and the inner region of V where V is positive. Such overlap 
persists dynamically and is responsible for producing a quadratic nonlinearity in the right 
hand side of the ^-equation. Similar observation can be made for the self-similar profiles 
for Initial Condition III, see Figured! 

As we can see from Figures [8] and [9J the rescaled profiles of the viscous nonlocal system 
is qualitatively similar to those of the inviscid nonlocal systems. This is to be expected since 
there is a logarithmic correction in the rescaling variable £ in the inviscid nonlocal system. 
Consequently, the viscous term can not dominate the nonlinear term in the nonlocal system. 
On the other hand, we observe that the profiles corresponding to the viscous nonlocal system 
are wider and more symmetric than those corresponding to the inviscid nonlocal system. 
This seems to make sense because the viscosity tends to smooth the singularity and make 
the profiles smoother and more symmetric. 

We have also performed a similar numerical study of the viscous nonlocal system with 
v = 0.01 for Initial Conditions II and III. We find that the viscous nonlocal system develops 
an asymptotically self-similar singularity in a finite time with the same blowup rate and 
self-similar scaling as the case of v = 0.001. 
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Figure 6: The inverse of ||u||oo (black) versus the asymptotic fit (red) for Initial Condition 
II with viscosity v = 0.001. The fitted blowup time is T = 0.833919962702315 and the 
scaling constant is C = 1.69630372479547. 




Time 



Figure 7: The inverse of ||u||oo (black) versus the asymptotic fit (red) for Initial Condition 
III with viscosity v = 0.001. The fitted blowup time is T = 0.617315651741129 and the 
scaling constant is C = 1.69150344092375. 
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Figure 8: Rescaled profiles U and V for Initial Condition II with viscosity v = 0.001 at 
t = 0.833917828434707, 0.83391976141767 and 0.833919943745501 respectively. The cor- 
responding maximum values of u are 794399, 8416207 and 89496701 respectively. Blue: 
profile of u; Red: profile of v. 
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Figure 9: Rescaled profiles U and V for Initial Condition III with viscosity v = 0.001 
at t = 0.617313605456105,0.617315459830455 and 0.617315633726175 respectively. The 
corresponding maximum values of u are 826395, 8808734 and 94072100 respectively. Blue: 
profile of u; Red: profile of v. 
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